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We show that the interplay between charging effects and the non-equilibrium spin accumulation 
has a dramatic effect in the current through an array of nanostructures attached to ferromagnetic 
electrodes. Large oscillations in the current as a function of bias voltage show up for parallel 
orientation of the electrodes' magnetizations. These oscillations originate in the inhomogeneity of 
the spin potentials through the array and correlate with oscillations in the spin accumulation. For 
antiparallel orientation the spin potential is homogeneous and the oscillations do not show up. This 
sensitivity results in a huge enhancement of the tunneling magnetoresistance as compared to the 
single-island case, and open new routes for improving the spintronic response of nanodevices. 

PACS numbers: 



A lot of effort has been devoted in the last two decades 
to understanding and controlling the interplay between 
magnetism and charge current because of promising ap- 
plications in electronic devices [IrE]- Interest is focused 
on systems with a large magnetoresistance (MR) , defined 
in terms of the resistance of the device for parallel (P) 
and antiparallel (AP) orientation of two ferromagnetic 
electrodes, MR = (Rap — Rp)/Rp. The tunneling mag- 
netoresistance (TMR) of magnetic tunnel junctions, early 
studied by Julliere[5], is controlled by the spin polariza- 
tion of the carriers p, assumed here to be equal in both 
electrodes: 
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In a magnetic double tunnel junction, with a metallic 
area inserted between the ferromagnetic electrodes, the 
TMR vanishes except if the spin flip processes are of lit- 
tle importance in the metal. Otherwise, the tunneling 
electrons lose memory of their initial spin. When the 
spin relaxation time is long, spin accumulation happens 
and induces a spin splitting of the chemical potential in 
the antiparallel configuration]!)]. The TMR is finite but 
only the half of Julliere's value. A further reduction of 
the magnetoresistance occurs with increasing number of 
metallic non-magnetic insertions, N, separated among 
themselves by tunnel barries: 
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where all the tunnel barriers separating the metallic in- 
sertions are supposed to be equivalent. The reduction of 
the TMR with the number of metallic regions is expected 
as it is only via the spin accumulation that the informa- 
tion about the relative magnetic orientation of the elec- 
trodes is transferred. This mechanism is indirect and less 
effective than the direct coupling of the two electrodes in 
a single tunnel junction]?]. 



If the dimensions of the metallic insertion are re- 
duced charging effects can no longer be neglected and 
new phenomena appear in the current [8] and in the 
magnetoresistance (9]. Due to the cost in energy to add an 
electron to the islands, at zero temperature, the current 
is blocked below a threshold voltage Vth- For a symmetri- 
cally biased array with nonmagnetic electrodes and short- 
range interactions Vth ~ 2Ai? c [TH] . Above this thresh- 
old the current is non-linear in voltage and frequently 
shows Coulomb staircase features. For a single island, a 
finite polarization of the electrodes can induce regions of 
negative differential conductance, as well as oscillations 
and changes of sign in the TMR, dependending on the 
specific set-up [T0HI5] . Charging effects are known to be 
enhanced in nanostructure arrays [8l IT6HT8] . however lit- 
tle is known about the spin transport in these systems. 
As the spin accumulation is the only way to transfer the 
spin information, a reduction of the TMR with increas- 
ing number of islands, similar to that in Eq.(2) could be 
expected, a priori. 

Here we show that, increasing the number of islands 
in between two ferromagnetic electrodes has a dramatic 
impact on the charge and spin transport. Large cur- 
rent oscillations show up in the I-V curves for parallel 
magnetic orientation, resulting in a strongly voltage de- 
pendent tunneling magnetoresistance. This behavior is 
a consequence of the interplay between charging effects 
and the spin accumulation through the array. The latter 
is homogenous along the array in the antiparallel con- 
figuration while it is inhomogenous and changes sign for 
parallel orientation. Opposite to what happens in the 
absence of charging effects, the TMR can be orders of 
magnitude larger than in the single particle case. 

We consider an array of N metallic nanostructures, in 
the following called islands, with charging energy E c and 
single particle level spacing 5 satisfying 8 -C K~bT -C E c . 
T is the temperature and Kg the Boltzmann constant. 
To allow for spin accumulation 8 is kept finite. The ar- 
ray is placed in between two ferromagnetic electrodes 
with spin polarization p. The islands are separated be- 
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tween themselves and from the ferromagnetic electrodes 
by equivalent tunnel junctions. The electronic interac- 
tions are assumed to be finite only when the charges 
are in the same conductor i.e. capacitive coupling be- 
tween different conductors vanishes, and charge disorder 
is absent. The electronic charge is taken equal to unity. 
Transport is treated at the sequential tunneling level with 
tunneling rates 

MABff) = ^^(A^TFl (3) 

Here AE a is the change in energy of electrons with spin 
a due to the tunneling process. The ferromagnetic po- 
larization of the electrodes enters via the tunneling re- 
sistance. R& is equal to 2Rt at the internal junctions 
separating two islands and to 2R T {1 ±p) _1 at the con- 
tact junctions between an electrode and the neighbor- 
ing island. Plus (minor) signs are assigned to majority 
(minority) spin carriers. The spin is conserved in the 
tunneling and in between tunneling events. Magneto- 
Coulomb [HlEl [20], spin-orbit effects, and the role of the 
metal insulator interface^!) in determining the spin po- 
larization of the carriers are neglected through all the 
paper. 

Under these assumptions the spin dependent potential 
at island i is <pi t<T = Ni i(7 5 + 2(A^ !(T + Ni^ a ) E c with i run- 
ning from 1 to N and N i tJ the number of excess electrons 
with spin a at island i. Correspondingly the spin accu- 
mulation at island i is given by A(f>i iCr — (iVj ;(7 — Ni^ a )5. 
The electrodes are maintained at spin-independent po- 
tentials 0o. a = V/2 and 4>N+i,a = —Vj2 and current 
flows perpendicular to the array axis. 

To calculate the current we use a Monte Carlo simula- 
tion which depends on the tunneling rates, as described 
in previous works [TBI US], At each iteration a single tun- 
neling event takes place. The time involved in this event 
depends on the tunneling rates of all the possible tunnel- 
ing processes. These rates are calculated at each itera- 
tion as they change when the spin and charge state of the 
system does. The values of the current given here corre- 
spond to stationary states. This computational method 
is appropriate for low temperatures when the electrons 
flow in the same direction in all the relevant tunneling 
processes. With increasing temperature the electrons be- 
come frequently stacked, oscillating back and forth be- 
tween two islands and computation is not possible. Cal- 
culations reported here correspond to p = 0.7, S = 10~ 5 
and KbT — 10~ 4 with energies given through all the 
paper in units of E c . The value of 5 does not affect the 
results as far as 5 -C KbT is satisfied. The non-zero tem- 
perature used to ensure 6 -C KbT allows a very small but 
finite current below threshold values. For simplicity, in 
most part of the discussion this finite temperature effect 
is neglected. 

Results for a single island are shown in Fig. 1. For par- 
allel arrangement the current can be described in terms 




FIG. 1: (a) Main figure: I-V curves of a single island placed 
in between two metallic (non-magnetic) or ferromagnetic elec- 
trodes for p = 0.7 with parallel and antiparallel orientation 
of the electrode magnetizations, see text. The corresponding 
TMR is shown as a function of voltage in the inset (a) where 
the Julliere's TMR is marked with a dashed line, (b) Main 
figure: Spin accumulation for a single island and antiparallel 
arrangement. Inset in (b): Sketch of the device and the spin 
accumulation for antiparallel orientation. The arrow in the 
island shows the spin with largest potential. 

of two independent spin channels and has the same value 
found with nonmagnetic electrodes, see Fig. la. For an- 
tiparallel orientation spin accumulation appears to equal- 
ize the ratios for entrance and exit from the island of spin 
up and down electrons, see inset in Fig. lb. Both the cur- 
rent and the spin accumulation increase monotonously 
with voltage above V t h as shown in Fig. 1. The inset in 
Fig. la shows the TMR corresponding to a single island 
as a function of the bias voltage V . When V ^> E c charg- 
ing effects are not important and the TMR equals half of 
Julliere's value, as expected for a double tunnel junction. 
As the voltage is reduced the TMR first decreases slightly 
and then increases until V t h where it suddenly jumps to 
zero. The maximum TMR, found at Vth equals Julliere's 
value given in Eq. 0. 

A very different behavior is found for N > 3, as the 
inner junctions between the islands can block the cur- 
rent flow. In the antiparallel configuration, as for non- 
magnetic electrodes, the current is strongly suppressed 
below V ~ Vth where a large jump is observed, see Figure 
2(a). Surprisingly for parallel orientation of the electrode 
magnetization the current show large oscillations below 
V t h ■ Within the range of voltages at which these oscilla- 
tions are observed, there are regions in which the current 
decreases with increasing V, corresponding to negative 
differential conductance. This negative conductance can 
be quite large, see the inset of Fig. 2 (a). As observed 
in Figure 2(b), the number of peaks observable in the 
current increases with the number of islands in the array. 

This unexpected behavior can be understood from an 
analysis of the potential drop through the array. In 




FIG. 2: (a) I-V curves for an N = 20 array with ferromag- 
netic electrodes in parallel and antiparallel arrangement. The 
I-V curve corresponding to nonmagnetic electrodes is plotted 
to serve as a reference. Inset: Differential conductance corre- 
sponding to N = 20 and parallel arrangement, (b) I-V curves 
for arrays of different sizes and parallel orientation of the elec- 
trode magnetizations. The number of peaks in the I-V curve 
increases with the size of the array. 



the case of metallic non-magnetic electrodes, the large 
threshold voltage originates in the lack of potential drop 
at the inner junctions which impede the flow of charges 
below Vth- This situation is typical of short-range in- 
teractions, as the voltage drops mostly at the junctions 
adjacent to source and drain, see discussion in jl8j. With 
increasing voltage charges are allowed to enter the array 
from the electrodes and create a charge gradient which 
provides the potential drops at the inner junctions nec- 
essary to allow the tunneling processes. Once a single 
charge occupies the first island, to increase the number 
of charges on it a voltage rise equal to 4E C is needed. Vth 
is the voltage necessary to put around N/2 charges onto 
the first island (the exact number depends on N being 
even or odd). 

Because of the spin accumulation when the electrodes 
are magnetic the potential profile is spin dependent. In 
the antiparallel configuration the spin-up and spin-down 
potentials change only very slightly when moving from 
an island to its neighboring ones, see Fig. 3(a). The spin 
accumulation is very homogeneous through the array. So, 
it barely contributes to the change in energy for tunneling 
at the inner junctions. As for non magnetic electrodes, 
this change in energy is controlled by the charge gradi- 
ent at the junctions. The position dependence of the spin 
accumulation is very different in the antiparallel and par- 
allel cases, see the sketches given at the top of Figure 3. 
For parallel orientation of the electrodes magnetization 
the spin accumulation changes sign as one moves from 
the first island to the end of the array, as seen in Fig. 
3(b), where we plot the spin potentials corresponding to 
N = 20 and V — 20E C as a function of the island po- 
sition. Spin down (up) potential increases (decreases) 
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FIG. 3: Top: Sketch of the spin accumulation in an N=6 ar- 
ray with antiparallel (left) and parallel configuration. Middle: 
Spin potentials as a function of island position at V = 20E C 
for an N — 20 array corresponding to antiparallel (a) and par- 
allel (b) orientation of the magnetization electrodes. Bottom: 
(c) and (d) show the voltage dependent spin polarization of 
the current and the spin accumulation in the first island for 
the same N = 20 array and parallel orientation. 



producing potential drops at the junctions between the 
islands which oppose (favor) the current flow. Because 
of this potential drop created by the spin accumulation 
the charge gradient needed to allow the flow of spin up 
electrons decreases. With the reduction of the number 
of charges which have to be accumulated at the first or 
last island to allow current, the threshold voltage is re- 
duced. Current is observed at smaller voltages for par- 
allel arrangement of the electrode magnetizations than 
in the antiparallel or nonmagnetic cases. With the po- 
tential drop at the inner junctions being spin dependent, 
naively one could expect a complete spin polarization of 
the current. As observed in Fig. 3(c), this does not hap- 
pen. The current polarization is not complete because 
the flow of spin-down electrons is allowed via the charge 
potential left by the spin-up ones. 

If, as discussed above, the current flows thanks to the 
spin accumulation, one could expect a decrease in the 
latter correlated with the regions of negative differential 
conductance. This is confirmed in Fig. 3(d) which shows 
oscillations in the spin accumulation at the same voltages 
as in the current. The decrease in the spin accumulation 
is due to the opening of a new transport channel only 
for spin down electrons. This is possible because the 
voltage drop necessary to permit the entrance of spin- 
down charges at the first island is smaller than the one 
corresponding to spin-up ones. At some V, being the 
first island occupied with n charges, spin-down electrons 
would be allowed to enter, while spin-up electrons will 
not. As spin accumulation happens to equilibrate the 
ratio between entrance and exit of spin up and down 
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FIG. 4: Tunneling magnetoresistance as a function of bias 
voltage for different array sizes. 



electrons, the opening of a new conduction channel for 
down electrons will reduce the spin accumulation. This 
is the origin of the decrease of current with increasing 
voltage, which lasts until this conduction channel is also 
opened for spin-up electrons. 

The differences in the current as a function of a mag- 
netic orientation have a dramatic effect on the TMR, 
plotted in Figure 4, which is strongly voltage dependent. 
Like the current it shows peaks in the parallel configura- 
tion. But the most impressive result are the extremely 
large values of the TMR observed. Even more than two 
orders of magnitude larger than in the single island case 
and than the value predicted by Eq. (2) which is recov- 
ered at large voltages (not shown). The largest values of 
the magnetoresistance depend on voltage and are largest 
for N = 3. This shows that there is a lot of room to 
improve the TMR response of electronic devices. 

In summary, we have studied the interplay between 
magnetism and charging effects in the transport through 
nanostructure arrays placed between ferromagnetic elec- 
trodes. We have found very non-linear I-V curves for 
parallel orientation of source and drain magnetizations, 
which include peaks in the current at voltages smaller 
than the metallic threshold voltage. This unusual de- 
pendence originates in the inhomogeneity of the spin po- 
tential through the array. The oscillations of the current 
correlate with oscillations of the TMR and of the spin 
accumulation which changes sign through the array. The 
TMR can reach values two orders of magnitude larger 
than those found in the single island case. Our results 
show that in the presence of charging effects a proper 
patterning can largely enhance the value of the magne- 
toresistance and will stimulate further experimental and 



theoretical studies. 
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